\(~\)
\(~\)

library(Seurat)
so = Read10X(data.dir = "/mnt/market/anclab-rstudio-server/home/swil0005/projects/Wnt11/raw_data/cellranger_mRNAcounts_asgz/")
so = CreateSeuratObject(counts = so , project = "Wnt11KO" , min.cells = 20)
so
## An object of class Seurat 
## 17039 features across 20586 samples within 1 assay 
## Active assay: RNA (17039 features, 0 variable features)

\(~\)
\(~\)

Cell Filtering

so[["percent.mt"]] <- PercentageFeatureSet(so, pattern = "^mt-")

quantile(so@meta.data$nCount_RNA , probs = seq(0, 1, by= 0.1))
##       0%      10%      20%      30%      40%      50%      60%      70% 
##   1047.0   1628.5   2651.0   3811.0   5708.0   7969.0  11349.0  14558.0 
##      80%      90%     100% 
##  18653.0  26267.5 200631.0
quantile(so@meta.data$nFeature_RNA , probs = seq(0, 1, by= 0.1))
##     0%    10%    20%    30%    40%    50%    60%    70%    80%    90%   100% 
##  242.0  631.0  704.0  889.0 1417.0 2326.0 3190.0 3670.0 4180.0 4954.5 9697.0
quantile(so@meta.data$percent.mt , probs = seq(0, 1, by= 0.1))
##          0%         10%         20%         30%         40%         50% 
##  0.04941627  0.40198602  0.81072654  1.61021127  1.95670073  2.19657747 
##         60%         70%         80%         90%        100% 
##  2.43974133  2.79772081  4.05519233 12.10581643 38.78822197
so <- subset(so, subset = nCount_RNA > 1500 & nCount_RNA < 20000 & nFeature_RNA > 500 & nFeature_RNA < 5000 & percent.mt < 5)
so[["RNA"]]
## Assay data with 17039 features for 11686 cells
## First 10 features:
##  Xkr4, Sox17, Mrpl15, Lypla1, Tcea1, Atp6v1h, Rb1cc1, 4732440D04Rik,
## Alkal1, Pcmtd1

\(~\)
\(~\)

Normalization

so <- NormalizeData(so , normalization.method = "LogNormalize" , scale.factor = 10000)

\(~\)
\(~\)

Finding the most variable genes

so <- FindVariableFeatures(so, selection.method = "vst", nfeatures = 2000)


LabelPoints(plot = VariableFeaturePlot(so), points = head(VariableFeatures(so), 10), repel = TRUE)

\(~\)
\(~\)

Mapping human gene symbols to mouse gene symbols

library(biomaRt)

s.genes = cc.genes$s.genes
g2m.genes = cc.genes$g2m.genes

# https://rjbioinformatics.com/2016/10/14/converting-mouse-to-human-gene-names-with-biomart-package/

human = useMart("ensembl", dataset = "hsapiens_gene_ensembl")
mouse = useMart("ensembl", dataset = "mmusculus_gene_ensembl")


s.genes = getLDS(attributes = c("hgnc_symbol"), filters = "hgnc_symbol", 
                 values = s.genes , mart = human,
                 attributesL = c("mgi_symbol"),martL = mouse, uniqueRows=T)
s.genes = s.genes$MGI.symbol


g2m.genes = getLDS(attributes = c("hgnc_symbol"), filters = "hgnc_symbol", 
                 values = g2m.genes , mart = human,
                 attributesL = c("mgi_symbol"),martL = mouse, uniqueRows=T)
g2m.genes = g2m.genes$MGI.symbol


# mgi_symbol: gene symboles for mouse
# hgnc_symbol: gene symbols for human

\(~\)
\(~\)

Cell Cycle scoring and Regression

# https://satijalab.org/seurat/archive/v3.2/cell_cycle_vignette.html

so <- ScaleData(so, features = rownames(so))

so = CellCycleScoring(so, s.features = s.genes, g2m.features = g2m.genes, set.ident = TRUE)
head(so[[]])
##                    orig.ident nCount_RNA nFeature_RNA percent.mt     S.Score
## AAACCCAAGGTCGTGA-1    Wnt11KO       7798          726  0.3334188 -0.02714792
## AAACCCACAACAAGAT-1    Wnt11KO      11151         3333  2.4123397 -0.12420683
## AAACCCACAGCTCGGT-1    Wnt11KO       1572          881  1.9083969  0.12395335
## AAACCCACATGAAGCG-1    Wnt11KO       4998          683  0.4201681 -0.04567019
## AAACCCAGTGTCCGTG-1    Wnt11KO      16636         3618  2.1399375  0.05585819
## AAACCCAGTGTGTACT-1    Wnt11KO       5014          862  0.3191065  0.04130828
##                       G2M.Score Phase old.ident
## AAACCCAAGGTCGTGA-1  0.003610455   G2M   Wnt11KO
## AAACCCACAACAAGAT-1 -0.244799795    G1   Wnt11KO
## AAACCCACAGCTCGGT-1  0.280071884   G2M   Wnt11KO
## AAACCCACATGAAGCG-1 -0.013205415    G1   Wnt11KO
## AAACCCAGTGTCCGTG-1  0.484910229   G2M   Wnt11KO
## AAACCCAGTGTGTACT-1  0.020404714     S   Wnt11KO
RidgePlot(so, features = c("Pcna", "Top2a", "Mcm6", "Mki67"), ncol = 2)

so <- RunPCA(so, features = VariableFeatures(so), nfeatures.print = 50)
DimPlot(so)

# Cells separate by cell-cycle phase
so <- RunPCA(so, features = c(s.genes, g2m.genes) , nfeatures.print = 10)
DimPlot(so)

# so <- ScaleData(so, vars.to.regress = c("S.Score", "G2M.Score"), features = rownames(so))
so = readRDS("Wnt11KO.rds")

# Now, a PCA on the variable genes no longer returns components associated with cell cycle
so <- RunPCA(so, features = VariableFeatures(so), nfeatures.print = 50)
DimPlot(so)

# When running a PCA on only cell cycle genes, cells no longer separate by cell-cycle phase
so <- RunPCA(so, features = c(s.genes, g2m.genes))
DimPlot(so)

\(~\)
\(~\)
\(~\)
\(~\)

Determining Dimentionality of Dataset

so = RunPCA(object = so , assay = "RNA" , features = NULL , npcs = 20 , rev.pca = FALSE ,
            verbose = TRUE, ndims.print = 1:3, nfeatures.print = 10, reduction.name = "pca", 
            reduction.key = "PC_", seed.use = 42, approx = T)



for(i in 1:20){
  DimHeatmap(so, dims = i, nfeatures = 40 , cells = 500 , reduction = "pca" , 
             disp.min = -6 , disp.max = 6 , balanced = T , projected = F, ncol = NULL,
             fast = TRUE, raster = TRUE, slot = "scale.data", assays = "RNA", combine = T)
}

# Selected PCs for clustering
# PCs: 1,2,3,4,5,6,7,8,13


ElbowPlot(so, ndims = 20, reduction = "pca")

\(~\)
\(~\)

Clustering with with embeddings

so1 <- FindNeighbors(so, return.neighbor = FALSE, compute.SNN = T,
                     prune.SNN = 1/15, nn.method = "annoy", n.trees = 200, annoy.metric = "euclidean", 
                     dims = c(1,2,3,4,5,6,7,8,13) , reduction = "pca" , verbose = F , do.plot = F )


# Cluster Determination
so1 = FindClusters(so1, modularity.fxn = 1 ,
                   initial.membership = NULL, node.sizes = NULL,
                   resolution = 1.5 , algorithm = 1,
                   n.start = 100, n.iter = 100, random.seed = 0,
                   group.singletons = T , temp.file.location = NULL,
                   edge.file.name = NULL, verbose = F)       



DimPlot(so1, reduction = "pca", label = T , dims = 1:2)

so1 = RunUMAP(so1, dims = 1:20)

head(so1[[]])
##                    orig.ident nCount_RNA nFeature_RNA percent.mt     S.Score
## AAACCCAAGGTCGTGA-1    Wnt11KO       7798          726  0.3334188 -0.03079170
## AAACCCACAACAAGAT-1    Wnt11KO      11151         3333  2.4123397 -0.11975109
## AAACCCACAGCTCGGT-1    Wnt11KO       1572          881  1.9083969  0.11639918
## AAACCCACATGAAGCG-1    Wnt11KO       4998          683  0.4201681 -0.04490011
## AAACCCAGTGTCCGTG-1    Wnt11KO      16636         3618  2.1399375  0.04813477
## AAACCCAGTGTGTACT-1    Wnt11KO       5014          862  0.3191065  0.03657714
##                       G2M.Score Phase old.ident RNA_snn_res.1.5 seurat_clusters
## AAACCCAAGGTCGTGA-1 -0.008638847    G1   Wnt11KO               2               2
## AAACCCACAACAAGAT-1 -0.257422641    G1   Wnt11KO               7               7
## AAACCCACAGCTCGGT-1  0.263503318   G2M   Wnt11KO               8               8
## AAACCCACATGAAGCG-1 -0.040737166    G1   Wnt11KO               1               1
## AAACCCAGTGTCCGTG-1  0.472502356   G2M   Wnt11KO              14              14
## AAACCCAGTGTGTACT-1  0.005532894     S   Wnt11KO               1               1
names(so1)
## [1] "RNA"     "RNA_nn"  "RNA_snn" "pca"     "umap"
DimPlot(so1, reduction = "umap", label = T)

\(~\)
\(~\)

Assigning cell type identity to clusters

features = "Six2"

VlnPlot(so1, features = features) + NoLegend()

RidgePlot(object = so1 , features = features) + NoLegend()

DotPlot(so1, features = features) + RotatedAxis()

FeaturePlot(so1, features = features)

\(~\)
\(~\)
\(~\)

Annotating Cluster 24:

Wnt4, C1qtnf12 (Fam132a), and Pax8 are the markers for Committing NP (N_Com).

S100g, Mal, Wfdc2, Aqp3, and Aqp2 are the markers for Distal nephron (N_Dist).

Wnt9b and Aqp2 are the markers for UE cortical collecting duct (UE_Ccd).

Upk3a, Upk1a, and Upk1b are the markers for UE Medullary collecting duct and Urothelium (UE_Mcd-Uro).

features = "Wnt4"

VlnPlot(so1, features = features) + NoLegend()

DotPlot(so1, features = features) + RotatedAxis()

FeaturePlot(so1, features = features)     

features = "S100g"

VlnPlot(so1, features = features) + NoLegend()

DotPlot(so1, features = features) + RotatedAxis()

FeaturePlot(so1, features = features)   

features = "Wnt9b"

VlnPlot(so1, features = features) + NoLegend()

DotPlot(so1, features = features) + RotatedAxis()

FeaturePlot(so1, features = features)      

features = "Upk3a"

VlnPlot(so1, features = features) + NoLegend()

DotPlot(so1, features = features) + RotatedAxis()

FeaturePlot(so1, features = features) 

\(~\)
\(~\)

### Nephron

# Nephron progenitor (N_NP), (Cited1, Crym, Six2 in clusters 4, 19)

# Committing NP/Committing nephron (PTA/RV/SSB)(N_Com), (Wnt4, in clusters 17,24), (C1qtnf12, inclusters 19,24,25,26), (Lef1 in cluster 17), (Pax8 in 24,26)

# Distal nephron/connecting segment (N_Dist), (S100g in 13,24,26), (Mal,Wfdc2,Aqp3,Aqp2 in 13,24), (Slc12a1 in 13)

# Nephron medial/Loop of Henle (N_LOH), No cluster.

# marker of Nephron proximal    (N_Prox), (Hnf4a, Lrp2, Aldob in 26)

# Nephron podocyte  (N_Pod), (Nphs1, Nphs2, Mafb, Synpo in 18)




### Ureteric Epithelium (Calb1 is the marker upregulated in 13)

# UE tip    (UE_Tip), (Ret,Wnt11,Crlf1 in 13). 

# UE cortical collecting duct   (UE_Ccd), (Wnt9b,Aqp2 in 13, 24), (Slit2 in 13).

# UE Medullary collecting duct and Urothelium   (UE_Mcd-Uro), (Upk3a, Upk1a in 24), (Upk1b in 24,25).





### Stroma (Pdgfra, Pdgfrb in 3,6,7,9,11,12,15,17,21)

# Cortical stroma   (S_CS), (Ntn1 in 6,12,21), (Foxd1 in 6,11,12,15,21), 

# Collecting-duct associated (S_CD), (Lef1 in 9,17), (Alx1 in 9,17)

# Stroma surrounding ureter S_UR, (Tbx18 in 11),

# Generic stromal markers   S,  




### Endothelium 

# Endothelial/vascular (Endo), (Pecam1 in clusters 20,23), (Cdh5, Flt1, Kdr in 20).




### Immune  

# Immune cells/myleoid cells    (Imm),  (Lyz2 in 10,14 (main),16,22,23), 
# (C1qa in 10,14(a little),16,22), (Fcer1g in 10, 14, 16, 22, 23)

\(~\)
\(~\)

Current.cluster.ids <- c(0:26)

New.cluster.ids <- c("Blood", "Blood", "Blood", "Stroma", "N_NP", "Blood", "S_CS" , "Stroma" , "8" ,
                     "S_CD" , "Imm" , "S_UR" , "S_CS" , "N_Dist/UE_Tip/UE_Ccd" , "Imm" , "S_CS" , 
                     "Imm" , "S_CD" , "N_Pod" , "N_NP" , "Endo" , "S_CS" , "Imm" , "Imm/Endo" , 
                     "N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd" , "N_Com/UE_Mcd-Uro" , "N_Prox/N_Com/N_Dist")

names(New.cluster.ids) <- levels(so1)
so1 <- RenameIdents(so1, New.cluster.ids)
DimPlot(so1, reduction = "umap", label = TRUE, pt.size = 0.5 , label.size = 3) + NoLegend()

\(~\)

so1 = AddMetaData(so1 , Idents(so1) , "Cell_Type")

head(so1[[]] , 20)
##                    orig.ident nCount_RNA nFeature_RNA percent.mt     S.Score
## AAACCCAAGGTCGTGA-1    Wnt11KO       7798          726  0.3334188 -0.03079170
## AAACCCACAACAAGAT-1    Wnt11KO      11151         3333  2.4123397 -0.11975109
## AAACCCACAGCTCGGT-1    Wnt11KO       1572          881  1.9083969  0.11639918
## AAACCCACATGAAGCG-1    Wnt11KO       4998          683  0.4201681 -0.04490011
## AAACCCAGTGTCCGTG-1    Wnt11KO      16636         3618  2.1399375  0.04813477
## AAACCCAGTGTGTACT-1    Wnt11KO       5014          862  0.3191065  0.03657714
## AAACCCATCAACACGT-1    Wnt11KO      15992         4076  3.0077539 -0.20074218
## AAACCCATCCGCAAAT-1    Wnt11KO       7300          724  0.2739726  0.01221797
## AAACCCATCTATCGTT-1    Wnt11KO       7010          713  1.0271041  0.03923562
## AAACCCATCTTGGCTC-1    Wnt11KO       1696         1023  2.2405660  0.07165411
## AAACGAAAGACCTGGA-1    Wnt11KO      12338         3359  1.9289998 -0.07249333
## AAACGAAAGCAGATAT-1    Wnt11KO       1721          806  1.1621150 -0.03952054
## AAACGAAAGCCTGACC-1    Wnt11KO       1695          744  1.2979351  0.19732527
## AAACGAAAGGAACGTC-1    Wnt11KO       2101          689  0.4759638  0.15713407
## AAACGAAAGGGCGAAG-1    Wnt11KO       7994         2793  2.4018014 -0.15401931
## AAACGAAAGGTCGTGA-1    Wnt11KO      12729         2934  2.3803912 -0.03005646
## AAACGAAAGTCATGAA-1    Wnt11KO      16758         3936  2.1780642 -0.13775837
## AAACGAACAAACTAAG-1    Wnt11KO      15533         3460  1.5386596  0.25344568
## AAACGAACAATGTCTG-1    Wnt11KO       3199         1473  1.1878712  0.20664185
## AAACGAAGTCCAGAAG-1    Wnt11KO      17039         3391  2.0775867 -0.19162184
##                       G2M.Score Phase old.ident RNA_snn_res.1.5 seurat_clusters
## AAACCCAAGGTCGTGA-1 -0.008638847    G1   Wnt11KO               2               2
## AAACCCACAACAAGAT-1 -0.257422641    G1   Wnt11KO               7               7
## AAACCCACAGCTCGGT-1  0.263503318   G2M   Wnt11KO               8               8
## AAACCCACATGAAGCG-1 -0.040737166    G1   Wnt11KO               1               1
## AAACCCAGTGTCCGTG-1  0.472502356   G2M   Wnt11KO              14              14
## AAACCCAGTGTGTACT-1  0.005532894     S   Wnt11KO               1               1
## AAACCCATCAACACGT-1 -0.281983255    G1   Wnt11KO              17              17
## AAACCCATCCGCAAAT-1 -0.096170188     S   Wnt11KO               2               2
## AAACCCATCTATCGTT-1 -0.046707833     S   Wnt11KO               0               0
## AAACCCATCTTGGCTC-1 -0.244101503     S   Wnt11KO               8               8
## AAACGAAAGACCTGGA-1  0.139565277   G2M   Wnt11KO              12              12
## AAACGAAAGCAGATAT-1  0.016302728   G2M   Wnt11KO               1               1
## AAACGAAAGCCTGACC-1  0.038377154     S   Wnt11KO               8               8
## AAACGAAAGGAACGTC-1 -0.010733195     S   Wnt11KO               1               1
## AAACGAAAGGGCGAAG-1 -0.313363478    G1   Wnt11KO              21              21
## AAACGAAAGGTCGTGA-1  0.432119583   G2M   Wnt11KO              20              20
## AAACGAAAGTCATGAA-1 -0.194367466    G1   Wnt11KO              13              13
## AAACGAACAAACTAAG-1 -0.191072581     S   Wnt11KO              16              16
## AAACGAACAATGTCTG-1 -0.074074612     S   Wnt11KO               8               8
## AAACGAAGTCCAGAAG-1 -0.295034457    G1   Wnt11KO               8               8
##                               Cell_Type
## AAACCCAAGGTCGTGA-1                Blood
## AAACCCACAACAAGAT-1               Stroma
## AAACCCACAGCTCGGT-1                    8
## AAACCCACATGAAGCG-1                Blood
## AAACCCAGTGTCCGTG-1                  Imm
## AAACCCAGTGTGTACT-1                Blood
## AAACCCATCAACACGT-1                 S_CD
## AAACCCATCCGCAAAT-1                Blood
## AAACCCATCTATCGTT-1                Blood
## AAACCCATCTTGGCTC-1                    8
## AAACGAAAGACCTGGA-1                 S_CS
## AAACGAAAGCAGATAT-1                Blood
## AAACGAAAGCCTGACC-1                    8
## AAACGAAAGGAACGTC-1                Blood
## AAACGAAAGGGCGAAG-1                 S_CS
## AAACGAAAGGTCGTGA-1                 Endo
## AAACGAAAGTCATGAA-1 N_Dist/UE_Tip/UE_Ccd
## AAACGAACAAACTAAG-1                  Imm
## AAACGAACAATGTCTG-1                    8
## AAACGAAGTCCAGAAG-1                    8

\(~\)
\(~\)

Finding differentially expressed features (DEFs) (cluster biomarkers).

so.markers <- FindAllMarkers(so1, slot = "data", min.pct = 0.1, verbose = F , 
                             logfc.threshold = 0.5 , test.use = "t", assay = "RNA", 
                             min.cells.feature = 10 , min.cells.group = 10)

# All Markers with absolute average log fold change larger than two
subset(so.markers , abs(so.markers$avg_log2FC) > 2 & so.markers$p_val_adj < 0.05)[,c(2,5,6)]
##                 avg_log2FC     p_val_adj                        cluster
## Hba-a1            5.009865  0.000000e+00                          Blood
## Hba-a2            5.001911  0.000000e+00                          Blood
## Hbb-bt            4.938371  0.000000e+00                          Blood
## Hbb-bs            4.919223  0.000000e+00                          Blood
## Alas2             3.486368  0.000000e+00                          Blood
## Mkrn1             2.993951  0.000000e+00                          Blood
## Bpgm              2.690968  0.000000e+00                          Blood
## Snca              2.677801  0.000000e+00                          Blood
## Fech              2.190977  0.000000e+00                          Blood
## Fam220a           2.070838  0.000000e+00                          Blood
## Slc25a37          2.019208  0.000000e+00                          Blood
## Cfl1             -2.042088  0.000000e+00                          Blood
## Col3a1           -2.048892  0.000000e+00                          Blood
## Lgals1           -2.090133  0.000000e+00                          Blood
## Col1a2           -2.109223  0.000000e+00                          Blood
## Meg3             -2.155789  0.000000e+00                          Blood
## Vim              -2.172572  0.000000e+00                          Blood
## Sparc            -2.204993  0.000000e+00                          Blood
## Cst3             -2.211365  0.000000e+00                          Blood
## Itm2b            -2.216454  0.000000e+00                          Blood
## Actb             -2.217154  0.000000e+00                          Blood
## Dcn              -2.246423  0.000000e+00                          Blood
## Malat1           -2.263107  0.000000e+00                          Blood
## mt-Cytb          -2.264501  0.000000e+00                          Blood
## Actg1            -2.291778  0.000000e+00                          Blood
## Col1a1           -2.350176  0.000000e+00                          Blood
## Cdkn1c           -2.438399  0.000000e+00                          Blood
## Tmsb4x           -2.997531  0.000000e+00                          Blood
## Apoe             -3.332003 4.147493e-249                          Blood
## Fcer1g           -2.276320 2.509331e-195                          Blood
## Tyrobp           -2.213460 8.715077e-191                          Blood
## Hbb-y             4.038364 1.749192e-137                          Blood
## Lyz2             -2.411889 1.689527e-133                          Blood
## C1qb             -2.222077 6.594791e-130                          Blood
## Hba-x             3.539685  4.496074e-62                          Blood
## Mkrn1.1          -2.232286  0.000000e+00                         Stroma
## Alas2.1          -2.541719  0.000000e+00                         Stroma
## Hbb-bt.1         -4.214194  0.000000e+00                         Stroma
## Hbb-bs.1         -4.319514  0.000000e+00                         Stroma
## Hba-a2.1         -4.396159  0.000000e+00                         Stroma
## Hba-a1.1         -4.424480  0.000000e+00                         Stroma
## Fcer1g.1         -2.003609 2.477138e-166                         Stroma
## Lyz2.1           -2.176337 2.058066e-146                         Stroma
## Mgp.1             2.009866 2.396092e-121                         Stroma
## Hbb-y.1          -4.876611  9.283181e-85                         Stroma
## Apoe.1           -2.749975  9.319485e-55                         Stroma
## Hba-x.1          -3.763994  9.329004e-50                         Stroma
## Fxyd6.1           2.077086  0.000000e+00                           N_NP
## Col1a2.2         -2.016121  0.000000e+00                           N_NP
## Col3a1.2         -2.087078  0.000000e+00                           N_NP
## Col1a1.2         -2.210267  0.000000e+00                           N_NP
## Alas2.2          -2.618031  0.000000e+00                           N_NP
## Tmsb4x.1         -3.172623  0.000000e+00                           N_NP
## Hbb-bt.2         -4.557309  0.000000e+00                           N_NP
## Hbb-bs.2         -4.608028  0.000000e+00                           N_NP
## Hba-a2.2         -4.699828  0.000000e+00                           N_NP
## Hba-a1.2         -4.783347  0.000000e+00                           N_NP
## Cited1.2          2.723417 1.683677e-319                           N_NP
## Mkrn1.2          -2.099821 4.108987e-279                           N_NP
## Apoe.2           -3.369404 3.731432e-243                           N_NP
## Fcer1g.2         -2.043260 1.830185e-202                           N_NP
## Dcn.2            -2.064448 1.588360e-199                           N_NP
## Lyz2.2           -2.176970 3.892392e-171                           N_NP
## C1qb.2           -2.018274 1.843467e-153                           N_NP
## Hbb-y.2          -4.665824  4.162694e-94                           N_NP
## Hba-x.2          -3.730137  5.668646e-41                           N_NP
## Foxd1.1           2.210584  0.000000e+00                           S_CS
## Aldh1a2.2         2.079697  0.000000e+00                           S_CS
## Mkrn1.3          -2.214426  0.000000e+00                           S_CS
## Alas2.3          -2.540206  0.000000e+00                           S_CS
## Hbb-bs.3         -4.185937  0.000000e+00                           S_CS
## Hbb-bt.3         -4.197440  0.000000e+00                           S_CS
## Hba-a2.3         -4.351547  0.000000e+00                           S_CS
## Hba-a1.3         -4.356825  0.000000e+00                           S_CS
## Meg3.2            2.070969 2.825017e-251                           S_CS
## Dlk1.2            2.285813 1.035762e-245                           S_CS
## Tyrobp.3         -2.030311 2.367143e-226                           S_CS
## Fcer1g.3         -2.078628 8.674893e-217                           S_CS
## Lyz2.3           -2.192568 3.745981e-157                           S_CS
## C1qb.3           -2.013512 1.427166e-135                           S_CS
## Hbb-y.3          -3.257407  1.847239e-53                           S_CS
## Apoe.3           -2.795011  1.748306e-49                           S_CS
## Hba-x.3          -2.853523  2.040742e-23                           S_CS
## Cdkn1c.2          2.686545 4.726188e-293                           S_CD
## Alas2.5          -2.389891 1.377525e-245                           S_CD
## Hba-a1.4         -4.034875 5.933338e-218                           S_CD
## Hba-a2.4         -4.042353 1.083976e-217                           S_CD
## Hbb-bt.5         -3.904486 2.539240e-217                           S_CD
## Hbb-bs.5         -3.833715 2.343454e-216                           S_CD
## Penk.3            2.604979 8.623620e-178                           S_CD
## Fcer1g.5         -2.023176 6.655842e-177                           S_CD
## Lyz2.5           -2.120207 9.388315e-114                           S_CD
## Hbb-y.4          -3.490904  1.661583e-27                           S_CD
## Apoe.5           -2.685152  4.851766e-23                           S_CD
## Hba-x.4          -2.755749  4.275659e-11                           S_CD
## Apoe.6            5.745254  0.000000e+00                            Imm
## Lyz2.6            4.898655  0.000000e+00                            Imm
## C1qb.6            4.811275  0.000000e+00                            Imm
## Fcer1g.6          4.660813  0.000000e+00                            Imm
## Tyrobp.6          4.386300  0.000000e+00                            Imm
## Coro1a.6          3.587922  0.000000e+00                            Imm
## Cd52.6            3.573172  0.000000e+00                            Imm
## Ctss.6            3.500510  0.000000e+00                            Imm
## Ctsc.6            3.261860  0.000000e+00                            Imm
## Ctsb.6            3.247724  0.000000e+00                            Imm
## Tmsb4x.5          3.221111  0.000000e+00                            Imm
## Lst1.6            3.209396  0.000000e+00                            Imm
## Ms4a6c.6          3.206885  0.000000e+00                            Imm
## Cyba.3            3.084055  0.000000e+00                            Imm
## Selenop.3         3.039459  0.000000e+00                            Imm
## Arhgdib.6         3.018058  0.000000e+00                            Imm
## Aif1.6            2.919618  0.000000e+00                            Imm
## Psap.2            2.855854  0.000000e+00                            Imm
## Lgmn.5            2.823587  0.000000e+00                            Imm
## Ptpn18.5          2.812157  0.000000e+00                            Imm
## Fcgr3.6           2.734368  0.000000e+00                            Imm
## Spi1.5            2.692691  0.000000e+00                            Imm
## Ctsd.1            2.684440  0.000000e+00                            Imm
## Cst3.4            2.678812  0.000000e+00                            Imm
## Rac2.6            2.671791  0.000000e+00                            Imm
## Laptm5.5          2.609907  0.000000e+00                            Imm
## Ctsh.5            2.572160  0.000000e+00                            Imm
## Ly86.5            2.557617  0.000000e+00                            Imm
## B2m.3             2.541774  0.000000e+00                            Imm
## Sat1.5            2.537128  0.000000e+00                            Imm
## Csf1r.5           2.531573  0.000000e+00                            Imm
## Arpc1b.2          2.499485  0.000000e+00                            Imm
## Ms4a6b.5          2.477358  0.000000e+00                            Imm
## Itm2b.3           2.420265  0.000000e+00                            Imm
## Fyb.5             2.390871  0.000000e+00                            Imm
## Gmfg.5            2.345757  0.000000e+00                            Imm
## Cd68.2            2.331817  0.000000e+00                            Imm
## Ctsz.2            2.329625  0.000000e+00                            Imm
## Unc93b1.2         2.302544  0.000000e+00                            Imm
## H2-D1.5           2.289112  0.000000e+00                            Imm
## Ckb.3             2.282515  0.000000e+00                            Imm
## Grn.1             2.208881  0.000000e+00                            Imm
## Hexa.1            2.175718  0.000000e+00                            Imm
## Rgs10.2           2.164107  0.000000e+00                            Imm
## Pld4.1            2.160643  0.000000e+00                            Imm
## Sh3bgrl3.2        2.152989  0.000000e+00                            Imm
## Lamp1.1           2.125554  0.000000e+00                            Imm
## Cotl1.2           2.108484  0.000000e+00                            Imm
## Clta.1            2.095699  0.000000e+00                            Imm
## Clec4a2.1         2.081758  0.000000e+00                            Imm
## Ehd4.1            2.078501  0.000000e+00                            Imm
## Tspo.3            2.061047  0.000000e+00                            Imm
## Psmb8.1           2.034493  0.000000e+00                            Imm
## Mdk.4            -2.186830  0.000000e+00                            Imm
## Sparc.3          -2.251869  0.000000e+00                            Imm
## Cdkn1c.3         -2.697607  0.000000e+00                            Imm
## C1qc.6            4.270142 1.168472e-316                            Imm
## Col1a1.5         -2.105762 1.330768e-311                            Imm
## Cx3cr1.1          2.239852 2.777859e-306                            Imm
## Hba-a1.5         -3.602695 5.153050e-292                            Imm
## Hba-a2.5         -3.576468 3.593782e-286                            Imm
## Hbb-bs.6         -3.503388 3.790010e-284                            Imm
## Alas2.6          -2.255738 1.697110e-282                            Imm
## C1qa.6            4.337589 3.379972e-280                            Imm
## Hbb-bt.6         -3.384944 2.564061e-268                            Imm
## Alox5ap.6         2.438275 1.453811e-256                            Imm
## Slfn2.1           2.004201 6.142255e-255                            Imm
## Dcn.4            -2.023512 1.164122e-211                            Imm
## Ly6e.2            2.001744 5.593573e-209                            Imm
## Fos.4             2.541716 2.365344e-203                            Imm
## Fcrls.1           2.295141 8.592411e-190                            Imm
## Mt1.5             2.158338 4.360272e-183                            Imm
## Dab2.2            2.155805 2.597052e-180                            Imm
## Pf4.6             4.001026 2.977531e-172                            Imm
## Ifitm3.3          2.782791 1.549683e-171                            Imm
## Ms4a7.1           2.081041 3.561917e-163                            Imm
## Mrc1              2.224318 7.683972e-145                            Imm
## Lgals3.1          2.463998 1.186633e-125                            Imm
## Ifi27l2a.1        2.220170 3.375187e-124                            Imm
## Hmox1.3           2.457087 1.475244e-106                            Imm
## Cd74              2.227840  3.547063e-52                            Imm
## Hbb-y.5          -4.356507  7.748176e-48                            Imm
## Hba-x.5          -3.434211  2.457970e-32                            Imm
## Col1a1.6          2.963275  0.000000e+00                           S_UR
## Col1a2.6          2.604663  0.000000e+00                           S_UR
## Sparc.4           2.369315  0.000000e+00                           S_UR
## Col3a1.6          2.552808 2.900421e-279                           S_UR
## Dcn.5             3.195732 2.058517e-259                           S_UR
## Meg3.5            2.254437 7.982242e-175                           S_UR
## Alas2.7          -2.412515 4.183164e-168                           S_UR
## Akap12.6          2.245148 1.563906e-158                           S_UR
## Cxcl12.4          2.371043 5.286460e-147                           S_UR
## Lum.4             3.188742 1.066729e-145                           S_UR
## Hbb-bs.7         -3.954088 3.901525e-144                           S_UR
## Hba-a1.6         -4.069833 6.856096e-142                           S_UR
## Hbb-bt.7         -3.921396 2.382666e-140                           S_UR
## Dlk1.6            2.247238 8.160804e-140                           S_UR
## Hba-a2.6         -4.050937 1.924021e-139                           S_UR
## Sfrp2.3           2.312653 6.549013e-134                           S_UR
## Mkrn1.7          -2.098201 4.846172e-127                           S_UR
## Lyz2.7           -2.075025  1.825218e-77                           S_UR
## Hbb-y.6          -2.502790  6.691056e-12                           S_UR
## Hba-x.6          -2.829514  5.347443e-06                           S_UR
## Epcam.7           3.468537 3.946723e-312           N_Dist/UE_Tip/UE_Ccd
## Wfdc2.6           3.437953 1.338285e-270           N_Dist/UE_Tip/UE_Ccd
## Ldhb.2            2.238744 1.157932e-254           N_Dist/UE_Tip/UE_Ccd
## Cldn7.1           2.586154 5.531850e-247           N_Dist/UE_Tip/UE_Ccd
## Vim.6            -2.410637 2.872840e-222           N_Dist/UE_Tip/UE_Ccd
## Col1a1.7         -2.269590 7.278929e-183           N_Dist/UE_Tip/UE_Ccd
## Col3a1.7         -2.117079 1.399181e-182           N_Dist/UE_Tip/UE_Ccd
## Col1a2.7         -2.050271 8.968457e-175           N_Dist/UE_Tip/UE_Ccd
## Krt18.6           2.524894 1.911047e-173           N_Dist/UE_Tip/UE_Ccd
## Lgals1.6         -2.364860 1.860429e-168           N_Dist/UE_Tip/UE_Ccd
## Calb1.7           4.093941 7.014694e-165           N_Dist/UE_Tip/UE_Ccd
## Dcn.6            -2.162377 2.647358e-159           N_Dist/UE_Tip/UE_Ccd
## Alas2.8          -2.471052 1.317094e-158           N_Dist/UE_Tip/UE_Ccd
## Cldn3             2.374805 8.619318e-156           N_Dist/UE_Tip/UE_Ccd
## Krt8.3            2.192433 7.328108e-155           N_Dist/UE_Tip/UE_Ccd
## Hbb-bt.8         -4.551041 9.181954e-149           N_Dist/UE_Tip/UE_Ccd
## Hba-a1.7         -4.492709 1.992368e-147           N_Dist/UE_Tip/UE_Ccd
## Hbb-bs.8         -4.219864 3.184453e-147           N_Dist/UE_Tip/UE_Ccd
## Hba-a2.7         -4.490026 2.130870e-146           N_Dist/UE_Tip/UE_Ccd
## Rprm              2.748766 1.448208e-124           N_Dist/UE_Tip/UE_Ccd
## Cdkn1c.5         -2.297807  1.910830e-98           N_Dist/UE_Tip/UE_Ccd
## Lyz2.8           -2.066784  1.228262e-62           N_Dist/UE_Tip/UE_Ccd
## Apoe.8           -2.883635  4.703977e-45           N_Dist/UE_Tip/UE_Ccd
## Hbb-y.7          -3.422496  2.162567e-33           N_Dist/UE_Tip/UE_Ccd
## Hba-x.7          -2.642879  3.408679e-09           N_Dist/UE_Tip/UE_Ccd
## S100g             2.047456  2.851832e-08           N_Dist/UE_Tip/UE_Ccd
## Igfbp7.5          2.515136  5.162524e-11                          N_Pod
## Hba-a2.8         -2.107410  2.210355e-10                          N_Pod
## Hba-a1.8         -2.074964  3.090650e-10                          N_Pod
## Mafb.2            3.063028  3.590388e-10                          N_Pod
## Hbb-bs.9         -2.002717  1.159542e-09                          N_Pod
## Pcp4              2.024515  1.121237e-08                          N_Pod
## Rasl11a.1         2.138632  2.426344e-08                          N_Pod
## R3hdml            2.656770  3.874996e-06                          N_Pod
## Apoe.9           -2.388812  4.481074e-02                          N_Pod
## Egfl7             3.051555  8.751735e-65                           Endo
## Crip2.2           2.349887  3.234954e-56                           Endo
## Fabp5.2           2.753293  4.083615e-56                           Endo
## Fabp4             4.609578  7.462791e-54                           Endo
## Icam2             2.397084  1.298368e-52                           Endo
## Emcn              2.894633  3.426105e-52                           Endo
## Cldn5.1           2.632344  4.089893e-49                           Endo
## Plvap             3.200724  1.884233e-47                           Endo
## Ramp2.4           2.269281  6.870409e-47                           Endo
## Gng11.2           2.163268  1.007255e-43                           Endo
## Fam167b           2.115451  3.655401e-43                           Endo
## Lyz2.10          -2.007078  3.782979e-18                           Endo
## Hba-a1.9         -2.301117  3.889220e-18                           Endo
## Hba-a2.9         -2.273174  8.589599e-18                           Endo
## Hbb-bt.10        -2.367982  1.790171e-17                           Endo
## Hbb-bs.10        -2.128467  1.834480e-17                           Endo
## Apoe.10          -2.178600  2.360740e-04                           Endo
## Coro1a.11         2.731706  8.312548e-61                       Imm/Endo
## Arhgdib.9         2.809156  7.675932e-59                       Imm/Endo
## Rac2.11           3.068671  2.047729e-53                       Imm/Endo
## Sh3bgrl3.5        2.053310  1.830868e-49                       Imm/Endo
## Gmfg.9            2.396527  1.660922e-37                       Imm/Endo
## Srgn.4            3.996566  1.586487e-35                       Imm/Endo
## Tyrobp.11         3.063740  1.431506e-34                       Imm/Endo
## Fcer1g.11         2.031083  2.911497e-34                       Imm/Endo
## Laptm5.9          2.037861  2.497048e-31                       Imm/Endo
## Tspo.5            2.046180  7.119600e-31                       Imm/Endo
## Cd52.11           2.795578  2.078738e-30                       Imm/Endo
## H2-D1.8           2.372063  3.636345e-30                       Imm/Endo
## Lcp1.2            2.217571  1.629972e-28                       Imm/Endo
## Alas2.11         -2.261315  3.901658e-25                       Imm/Endo
## Fxyd5.2           3.127217  5.593747e-25                       Imm/Endo
## Man2b1.2          2.108062  1.040911e-19                       Imm/Endo
## Emb.2             2.007622  1.244493e-19                       Imm/Endo
## Hba-a2.10        -3.359457  1.382025e-19                       Imm/Endo
## Fcgr3.11          2.035771  3.357833e-19                       Imm/Endo
## Hba-a1.10        -3.370828  4.082851e-19                       Imm/Endo
## Hbb-bt.11        -3.378443  1.510421e-18                       Imm/Endo
## Ptprcap           2.320070  5.492539e-18                       Imm/Endo
## Hbb-bs.11        -3.134245  8.107288e-18                       Imm/Endo
## Alox5ap.11        3.542391  2.120266e-17                       Imm/Endo
## Hdc               3.564075  4.998452e-14                       Imm/Endo
## Ifitm1.5          6.049862  1.366749e-13                       Imm/Endo
## Cd9.1             2.033821  2.527257e-13                       Imm/Endo
## Ltb               2.184038  8.656467e-13                       Imm/Endo
## Plac8.6           2.310042  1.574734e-12                       Imm/Endo
## Lgals3.2          2.276612  4.565328e-12                       Imm/Endo
## Pglyrp1           2.968278  9.102910e-10                       Imm/Endo
## Cpa3              4.401511  2.154801e-09                       Imm/Endo
## Cyp11a1           2.382940  8.700792e-09                       Imm/Endo
## Cd69              2.116353  2.276980e-07                       Imm/Endo
## Cd200r3           2.171432  4.762573e-07                       Imm/Endo
## Ccl9.1            2.723936  1.054960e-06                       Imm/Endo
## Ccl3.1            2.599118  2.081037e-05                       Imm/Endo
## Mcpt8             3.376087  4.712338e-04                       Imm/Endo
## Krt18.8           4.209438  3.437371e-80 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Epcam.11          3.125237  1.942720e-78 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Krt8.4            4.045806  1.024845e-73 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Alas2.12         -2.661593  5.231427e-71 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Cldn7.3           3.280594  4.502295e-69 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Spint2.2          2.404951  9.420929e-66 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Wfdc2.9           3.015752  1.573198e-60 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Tacstd2.1         3.422996  4.103121e-53 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Bcam.2            2.469606  3.258920e-51 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Col1a2.11        -2.114803  4.988994e-50 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Krt7.2            3.319589  1.572614e-48 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Atp1b1.3          2.264765  3.020460e-48 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Hbb-bt.12        -5.099266  7.968804e-47 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Lgals1.10        -2.614273  1.016680e-44 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Vim.8            -2.371101  1.127828e-44 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Hba-a1.11        -4.945765  1.588618e-44 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Krt19             3.707714  6.416738e-44 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Hba-a2.11        -4.982656  1.626493e-43 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Mest.10          -2.184883  2.313260e-43 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Hbb-bs.12        -4.666126  2.535880e-43 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Tspan8            2.936064  3.965445e-43 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Cldn4.1           2.244381  4.109752e-43 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Col1a1.11        -2.267846  1.027212e-42 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Lyz2.12          -2.131684  8.184368e-42 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Foxq1             2.647323  3.888820e-41 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Col3a1.11        -2.082213  6.245223e-41 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Hspb1.1           2.579797  2.757308e-38 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Mal.1             2.477846  5.447142e-38 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Cystm1            2.077164  2.978615e-37 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Wnt4.1            2.175289  5.372146e-34 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Pdzk1ip1.1        2.933248  1.845267e-32 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Dcn.10           -2.097509  1.388950e-31 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Apoe.11          -3.355543  3.290511e-28 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Cldn3.1           2.177516  1.069880e-26 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Fxyd4             3.719893  4.311397e-25 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Cdh16.2           2.275180  9.254549e-25 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Fxyd3             3.074420  6.887786e-24 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Id1.6             2.145059  1.767947e-23 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Sfn               2.090866  2.120981e-20 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## S100a6.10         2.074633  1.483915e-16 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Cdkn1c.9         -2.041002  3.322340e-13 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Upk3a             2.534381  4.385484e-12 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Hbb-y.10         -5.007761  7.326094e-09 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Hba-x.9          -3.880112  1.676136e-07 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Sparc.8           3.346582  1.980169e-57               N_Com/UE_Mcd-Uro
## Col3a1.12         2.673101  2.064136e-54               N_Com/UE_Mcd-Uro
## Sfrp2.8           3.431200  4.023557e-51               N_Com/UE_Mcd-Uro
## Col1a1.12         2.735051  5.039510e-49               N_Com/UE_Mcd-Uro
## Dlk1.12           3.566384  8.619757e-46               N_Com/UE_Mcd-Uro
## Fbln2.2           3.366598  5.894116e-45               N_Com/UE_Mcd-Uro
## Cryab.1           2.946992  9.215705e-45               N_Com/UE_Mcd-Uro
## Krt19.1           3.330028  1.886138e-44               N_Com/UE_Mcd-Uro
## Upk3b             3.784039  4.428858e-42               N_Com/UE_Mcd-Uro
## Crip1.8           3.226252  5.677367e-40               N_Com/UE_Mcd-Uro
## Meg3.10           2.424259  6.181061e-39               N_Com/UE_Mcd-Uro
## Slc9a3r1.4        2.503056  6.862073e-38               N_Com/UE_Mcd-Uro
## Hspb1.2           3.140294  9.232364e-38               N_Com/UE_Mcd-Uro
## Lyz2.13          -2.137057  1.161133e-37               N_Com/UE_Mcd-Uro
## Col14a1.10        2.756124  1.945911e-37               N_Com/UE_Mcd-Uro
## Cavin3.4          2.155528  3.661515e-36               N_Com/UE_Mcd-Uro
## Col1a2.12         2.168213  4.391321e-36               N_Com/UE_Mcd-Uro
## Krt7.3            2.819083  1.586142e-35               N_Com/UE_Mcd-Uro
## S100a6.11         3.476857  2.269283e-34               N_Com/UE_Mcd-Uro
## Tmem108           2.200778  1.831667e-30               N_Com/UE_Mcd-Uro
## Itm2a.9           3.388282  1.741411e-29               N_Com/UE_Mcd-Uro
## Hbb-bt.13        -4.608916  3.334300e-28               N_Com/UE_Mcd-Uro
## Alas2.13         -2.393809  1.227427e-27               N_Com/UE_Mcd-Uro
## Rspo1.1           2.009231  1.479489e-27               N_Com/UE_Mcd-Uro
## Hbb-bs.13        -4.124014  4.626799e-27               N_Com/UE_Mcd-Uro
## Hba-a1.12        -4.389347  4.757441e-27               N_Com/UE_Mcd-Uro
## Hba-a2.12        -4.379492  6.346572e-27               N_Com/UE_Mcd-Uro
## Eln               2.327724  2.249836e-25               N_Com/UE_Mcd-Uro
## Rarres2.10        2.427102  6.673742e-23               N_Com/UE_Mcd-Uro
## Aebp1.3           2.045606  1.125318e-22               N_Com/UE_Mcd-Uro
## 1500015O10Rik.1   2.264760  2.178616e-22               N_Com/UE_Mcd-Uro
## Igfbp6.6          2.247473  2.708995e-12               N_Com/UE_Mcd-Uro
## Fmo2              3.033400  9.313511e-12               N_Com/UE_Mcd-Uro
## Hba-x.10         -3.960187  1.432245e-11               N_Com/UE_Mcd-Uro
## Hbb-y.11         -5.118240  5.472239e-10               N_Com/UE_Mcd-Uro
## Stmn2.2           3.272816  8.936980e-08               N_Com/UE_Mcd-Uro
## Cyba.5            2.354472  7.533231e-35            N_Prox/N_Com/N_Dist
## Fxyd2.1           5.291368  2.245865e-34            N_Prox/N_Com/N_Dist
## Ly6a.3            4.212719  5.478547e-33            N_Prox/N_Com/N_Dist
## Lrpap1.2          2.651403  1.108613e-31            N_Prox/N_Com/N_Dist
## Calml4.2          3.179583  9.982252e-31            N_Prox/N_Com/N_Dist
## Slc9a3r1.5        2.215087  1.068975e-30            N_Prox/N_Com/N_Dist
## Mpc2.4            2.486508  1.197684e-30            N_Prox/N_Com/N_Dist
## Pdzk1             2.746115  1.829981e-29            N_Prox/N_Com/N_Dist
## Rida.1            3.037432  1.697521e-28            N_Prox/N_Com/N_Dist
## Chchd10.4         2.554388  1.904516e-27            N_Prox/N_Com/N_Dist
## Amn               2.654958  1.551127e-26            N_Prox/N_Com/N_Dist
## Mdh1.4            2.714683  2.462354e-26            N_Prox/N_Com/N_Dist
## Igfbp7.10         2.365169  3.596736e-26            N_Prox/N_Com/N_Dist
## Folr1.2           2.530971  9.119426e-26            N_Prox/N_Com/N_Dist
## Ezr.4             2.021599  9.895440e-25            N_Prox/N_Com/N_Dist
## Atp1b1.5          2.674680  1.059714e-24            N_Prox/N_Com/N_Dist
## Akr7a5.4          2.170087  1.805664e-24            N_Prox/N_Com/N_Dist
## Npnt.4            2.169077  1.910573e-24            N_Prox/N_Com/N_Dist
## Slc39a5           2.070364  3.003691e-24            N_Prox/N_Com/N_Dist
## Ggt1              3.105118  3.576864e-24            N_Prox/N_Com/N_Dist
## Ldhb.7            2.332231  6.241133e-24            N_Prox/N_Com/N_Dist
## Pcbd1.3           2.043439  2.726802e-23            N_Prox/N_Com/N_Dist
## Cpn1.1            2.623152  3.193002e-23            N_Prox/N_Com/N_Dist
## Phgdh.8           2.884691  1.248214e-22            N_Prox/N_Com/N_Dist
## Lgals1.12        -2.742854  1.309490e-21            N_Prox/N_Com/N_Dist
## 0610005C13Rik     2.225278  1.584856e-21            N_Prox/N_Com/N_Dist
## Clu.4             2.073325  4.494623e-21            N_Prox/N_Com/N_Dist
## Alas2.14         -2.586118  6.791025e-21            N_Prox/N_Com/N_Dist
## Glud1.3           2.132607  2.785794e-20            N_Prox/N_Com/N_Dist
## Hba-a2.13        -5.321656  6.168150e-20            N_Prox/N_Com/N_Dist
## Ttc36             4.087390  2.159100e-19            N_Prox/N_Com/N_Dist
## Tmsb4x.8         -3.405344  2.436396e-19            N_Prox/N_Com/N_Dist
## Hbb-bs.14        -5.159188  3.769772e-19            N_Prox/N_Com/N_Dist
## Dcn.12           -2.169356  6.936446e-19            N_Prox/N_Com/N_Dist
## Cyb5a.5           2.061208  1.193316e-18            N_Prox/N_Com/N_Dist
## Col1a2.13        -2.067640  4.111726e-18            N_Prox/N_Com/N_Dist
## Aldob             3.945110  5.924123e-18            N_Prox/N_Com/N_Dist
## Slc34a1           3.340307  6.437418e-18            N_Prox/N_Com/N_Dist
## Hbb-bt.14        -5.166921  1.171618e-17            N_Prox/N_Com/N_Dist
## Hba-a1.13        -5.129147  3.504662e-17            N_Prox/N_Com/N_Dist
## Vim.10           -2.263959  2.836193e-16            N_Prox/N_Com/N_Dist
## Keg1              2.088222  2.979535e-16            N_Prox/N_Com/N_Dist
## Hmgcs1.2          2.009830  1.987327e-15            N_Prox/N_Com/N_Dist
## Spink1            2.402158  9.335577e-15            N_Prox/N_Com/N_Dist
## Sult1d1           3.400172  9.478032e-15            N_Prox/N_Com/N_Dist
## Mep1a             2.132565  1.937744e-14            N_Prox/N_Com/N_Dist
## Col3a1.13        -2.021109  2.335123e-14            N_Prox/N_Com/N_Dist
## Mest.12          -2.010056  3.614430e-14            N_Prox/N_Com/N_Dist
## Mkrn1.14         -2.139525  4.679515e-13            N_Prox/N_Com/N_Dist
## Tmem174           2.075951  2.329875e-12            N_Prox/N_Com/N_Dist
## Ass1.1            2.253869  4.666053e-12            N_Prox/N_Com/N_Dist
## Col1a1.13        -2.133876  5.454298e-12            N_Prox/N_Com/N_Dist
## Mt1.10            3.561380  2.923428e-10            N_Prox/N_Com/N_Dist
## Sostdc1.2         2.013489  9.531778e-10            N_Prox/N_Com/N_Dist
## Fabp3             2.024052  6.916508e-09            N_Prox/N_Com/N_Dist
## Cdkn1c.11        -2.265757  6.929272e-09            N_Prox/N_Com/N_Dist
## Hdc.1             2.094512  1.997294e-08            N_Prox/N_Com/N_Dist
## Mt2               2.434612  6.027492e-08            N_Prox/N_Com/N_Dist
## Spp2              2.217776  9.397472e-08            N_Prox/N_Com/N_Dist
## Hbb-y.12         -5.097419  1.151633e-04            N_Prox/N_Com/N_Dist

\(~\)
\(~\)
\(~\)
\(~\)